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Dynamics  of  Electron  Transfer  for  a  Nonsuperexchange 

Coherent  Mechanism.  I. 

R.  A.  Marcus  and  R.  Almeida^ 

Noyes  Laboratory  of  Chemical  Physics* 

California  Institute  of  Technology 
Pasadena,  California  91125 

Abstract 

In  addition  to  mechanisms  such  as  superexchange  and  a  chemical  intermediate 
mechanism  for  electron  transfer  from  a  donor  D*  to  an  acceptor  A  via  a  molecular 
bridge  B,  a  third  possibility  occurs  when  the  I  dectronic  coupling  is  very  strong 
and  the  D*B  and  D+B~  states  have  energies  moderately  close  to  each  other.  This 
mechanism  is  discussed  here.  Like  superexchange,  it  is  a  coherent  one,  in  contrast  to 
the  chemical  intermediate  mechanism,  where  the  transfer  is  sequential  and 
incoherent.  The  dynamics  of  the  mechanism  are  discussed,  particularly  the 
maximum  population  of  "B~”  and  the  question  of  whether  an  effective  rate  constant 
for  its  disappearance  can  be  considerably  larger  than  the  maximum  adiabatic  rate 
constant.  There  are,  as  yet,  no  experimental  data  on  the  mechanism,  though  the 
synthesis  of  suitable  D*BA’s  may  permit  its  observation. 

In  the  treatment  three  collective  nuclear  coordinates  are  introduced,  permitting 
independent  reorganization  energies  for  each  reactive  center.  With  certain 
approximations,  namely,  equal  vibration  frequencies  and  a  nonadiabatic  first  step, 
the  problem  is  reduced  analytically  to  a  one-coordinate  one,  which  can  be  readily 
treated  numerically.  One  rough  but  simple  analytical  result  for  the  latter  is  also 
given. 
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I.  Introduction 

Electron  transfers  from  a  donor  D*  to  an  acceptor  A  via  a  bridge  B  may  occur  by 
one  of  several  mechanisms.  The  most  common  one  is  that  of  superexchange,  in  which 
the  transferring  electron  makes  use  of  a  virtual  state,  an  orbital  quite  different  in 
energy  from  that  of  the  donor  D*  or  the  acceptor  A.  A  second  mechanism,  which  could 
occur  when  the  relevant  bridge  orbital  is  close  enough  in  energy  to  the  D*  one,  is  for 
the  electron  transfer  to  occur  via  a  chemical  intermediate  B",  whose  presence  could 
be  detected  experimentally.  The  kinetics  of  this  transfer  can  be  treated 
straightforwardly  using  two  consecutive  steps. 

A  third  possibility  is  considered  in  the  present  paper,  one  for  which  there  are  as 
yet  no  experiments,  but  for  which  suitable  D*BA  systems  might  be  synthesized.  In 
this  mechanism  the  relevant  B  orbital  is  readily  accessible  energetically,  as  in  the 
second  mechanism  above,  but  B  and  A  are  so  strongly  coupled  electronically  that  the 
entire  transfer  occurs  coherently,  rather  than  incoherently  in  two  successive  steps. 

One  question  which  arises  concerns  the  maximum  population  of  "B-"  in  this  third 
mechanism.  In  the  case  of  the  second  (sequential)  mechanism  this  maximum 
population  B-^  is  given  by1 


B"  = 

max 


(l.i) 


where  kx  and  k2  are  the  rate  constants  of  the  D*-»B  and  the  B--»A  electron  transfers, 
respectively.  The  question  which  arises  is  whether  or  not  for  the  third  mechanism  eq. 
(1.1)  can  still  approximately  represent  the  data  but  with  a  considerably  enhanced 
value  of  k2 ,  a  value  substantially  greater  than  the  maximum  k2  for  an  adiabatic 
electron  transfer.  (In  this  case  B-^  would  be  much  smaller  than  for  a  sequential 


mechanism.)  The  answer,  as  we  shall  see  in  a  subsequent  article  on  numerical 
results,  is  in  the  affirmative.  Presumably  by  a  suitable  choice  of  D*BA  systems,  e.g., 
a  D*  weakly  coupled  to  B,  a  suitable  choice  of  B  with  an  energy  level  close  to  or  far 
above  that  of  D*  and  a  suitable  choice  of  the  electronic  coupling  between  B  and  A 
(e.g.,  an  aromatic  or  other  appropriate  group  provides  strong  coupling,  a  nonbornyl 
group  relatively  weaker  coupling)  it  may  be  possible  to  construct  a  series  of  systems, 
different  members  of  which  proceed  by  one  or  the  other  of  the  three  mechanisms. 

In  the  present  paper  we  set  up  a  model  for  performing  calculations  of  B-max  and 
other  properties.  In  this  model  there  should  preferably,  for  the  purpose  cited,  be  at 
least  three  collective  coordinates,2  each  having  its  own  reorganization  parameter  and 
each  associated  with  an  equilibrium  vibrational  displacement  in  a  reactive  center,  D, 
B  and  A.  There  should  also  be,  in  a  quantum  mechanical  treatment  of  the  dynamics 
of  this  coherent  transfer,  some  mechanism  for  dissipation,  at  least  towards  the  end  of 
the  transfer  process,  so  that  a  wavepacket  describing  the  dynamics  is  not  spuriously 
reflected  back  and  forth.  In  practice,  this  dissipation  is  provided  by  the  many 
coordinates  in  the  actual  system,  which  permit  the  reaction  products  to  become 
thermalized  before  any  reverse  reaction  occurs.  (To  design  a  theory  which  in  one 
limit  would  yield  two  incoherently  connected  consecutive  steps,  a  dissipative  term 
accompanying  the  motion  in  the  D+B~  A  stage  would  also  be  included.) 

Rather  than  solve  the  Schrbdinger  equation  with  these  three  collective 
coordinates  purely  numerically,  although  that  would  be  useful,  we  introduce  in  the 
present  paper  some  approximations  such  that  the  coherent  problem  can  be  reduced 
analytically  from  three  coordinates  to  one.  Then  the  latter  can  then  either  be  treated 
numerically  (e.g.,  using  a  one-dimensional  Fast  Fourier  Transform3)  or,  in  a  more 
approximate  way,  analytically.  A  rather  rough  estimate  is  given  in  Section  V. 

A  coordinate  system  and  the  Hamiltonian  is  set  up  in  Section  II  and  an 
approximation  is  introduced  (all  three  vibrational  frequencies  equal)  which,  when 


followed  by  a  rotation  of  the  coordinates,  permits  a  reduction  to  a  two-coordinate 
problem.  In  section  HI  the  transfer  from  D*  to  the  BA  coupled  system  is  treated 
nonadiabatically  (Golden  Rule).  To  calculate  the  time  evolution  of  B",  B"(f),  the 
subsequent  dynamics  are  investigated  (section  IV)  in  which  the  nonadiabatic  step 
deposits  a  wavepacket  on  the  BA  potential  energy  surface.  This  wavepacket, 
continuously  replenished  from  the  decreasing  concentration  of  D*,  evolves  in  time  on 
the  BA  surface,  and  has  a  nuclear  motion  on  BA  which  is  seen  in  section  IV  to  be 
separable.  The  treatment  of  the  coherent  dynamics  has  thereby  been  reduced  to  a 
one-coordinate  problem.  The  energy  partitioning  and  its  effect  on  B "(f)  are  discussed 
in  section  V.  Numerical  results  are  given  in  a  subsequent  paper.4 

II.  Coordinates  and  Hamiltonian 

Diabatic  electronic  wavefunctions  xpi  are  introduced  to  describe  the  three 
electronic  configurations  D*BA,  D  +  B~A,  and  D  +  BA~,  corresponding  to  the 
transferring  electron  being  on  D(i  =  l),  B(i  =  2)  and  A(i  =  3).  We  let  4>f  denote  the 
corresponding  time-dependent  nuclear  wavefunctions,  to  be  determined  by  solution  of 
the  Schrodinger  equation.  The  electronic-nuclear  wavefunction  qr(x)  at  time  x  is  now 
given  by 

3 

V|»(x)  =  ^  (*)  Vi  ■  (2.1) 

t  =  l 

The  Hamiltonian  contains  a  nuclear  kinetic  energy  operator  K  and  several 
nuclear-electronic  terms,  Hel.  Introduction  of  ip(x)  into  the  Schrodinger  equation 
H^-ih  dq//dx,  multiplication  by  qr*  and  integration  over  the  electronic  coordinates 
yields  a  standard  result  for  the  nuclear  motion: 


-4- 


(2.3) 


where  Hy  denotes  the  matrix  element 

IV 

We  next  introduce  three  collective  vibrational  coordinates,  qv  qv  g3,  one 

collective  vibration  per  center.2  For  simplicity,  a  common  vibration  frequency  co/2n  is 

* 

assumed  for  these  vibrations  (m  =  (ki  /y^)  ,  where  kt  is  the  force  constant  and  p(  a 

* 

reduced  mass),  and  dimensionless  coordinates  qt  =  (ptM/A)  qx  are  introduced. 
Equation  (2.2)  now  becomes 


—  + 


y  h..$. 


(2.4) 


where  t  is  a  dimensionless  time  variable  cot  and  Hi}  denotes  the  Hi}  in  eq.  (2.2)  divided 
by  Am.  In  actual  fact,  these  three  vibrational  frequencies  m  will  not  be  equal,  but  we 
are  more  interested  here  in  obtaining  some  insight  into  what  magnitudes  to  expect 
for  B~U)  for  this  mechanism,  rather  than  in  an  accurate  numerical  solution. 

For  each  of  the  three  electronic  configurations  D*BA,  D+B~  A  and  D+BA“,  a  set  of 
three  equilibrium  values  of  the  q\ s  is  defined.  When  the  electron  is  on  D,  namely 
when  is  large,  the  equilibrium  values  of  (qt,  q2,  q3)  are  denoted  by  (cjj.O.O),  while 
when  the  electron  is  on  B,  they  are  taken  to  be  (0,a2,0)  and,  when  on  A,  (0,0,a3)  (Fig. 
1).  Thus,  with  this  choice  we  have,  for  example,  q{  =  ctj  at  the  equilibrium  position  for 
the  vibration  qx  in  the  D*  system,  and  qx  =0  at  the  equilibrium  position  for  qx  when 
the  D  is  in  the  form  D  +  .  Using  a  harmonic  approximation  the  diagonal  matrix 
elements  Hu  are  now 


(2.5) 


Hn  =  Kq^-aJ2  +  ±q2  +  * q 2 

H22  =  ft*  +  *<<72  “a  2?  +  *<?32  +  ^12  f 


#33  =*^t2  +  *<722  +  *<<?3-a3)2  +  ^#13 


where  the  AEU’ s  denote  the  (equilibrium)  energy  differences  of  the  first  and  the  tth 
diabatic  electronic  configurations,  divided  by  ha.  The  off-diagonal  HJ s  will  be  taken 
to  be  independent  of  the  q’s. 

To  solve  eqs.  (2.4)  -  (2.5)  we  first  introduce  a  new  Cartesian  coordinate  system  (x, 
y,  z),  which  permits  a  separation  of  the  variable  x  from  the  {y,  z)  pair.  The  new  axes 
are  obtained  by  a  rotation  of  the  (qv  q2,  q3 )  axes,  as  in  Fig.  1.  Customarily,  three 
rotations  are  used  for  rotating  the  axes,  as  in  the  definition  of  the  Euler  angles.5 
However,  for  our  purposes  it  suffices  to  introduce  only  two:  The  rotations  are  such 
that  the  final  x-axis  is  perpendicular  to  the  aja2a3  plane  and  the  final  z-axis  is 
parallel  to  the  a2a3  line. 

We  first  rotate  the  axes  through  an  angle  $  about  the  qt-axis,  such  that  the  new 
q3-coordinate  (z  in  Fig.  1)  is  parallel  to  the  line  a2a3,  and  then  rotate  the  axes  through 
an  angle  0  about  this  z-axis,  such  that  the  new  qj-axis  (x)  is  perpendicular  to  the 
plane  ata2a3.  We  thus  have 


M 

^  COS0 

sin  0 

y 

= 

—  sin  0 

cos  0 

1 

\° 

0 

0 

cos$ 
—  sin  <J> 


where  in  Fig.  1  cos  $  =  0P/a2,  sin  $  =  0 P/a3,  cos  0  =  OQ/aj  and  sin  0  =  OQ/OP.  From  the 
latter  results  we  have  OP = A2  and  OQ  =  A  t ,  where 
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(2.7) 


H22(^)  =  /f33(>)  =  iCy-A1A2/a1)2 
H22(z)  =  -Hz  +  a^Aj/ag)2  +  &En 
H33(z)  =  Uz-  a3'A2/a2)2  +  A£,  3 
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With  the  change  of  coordinates  from  (qv  q2,  q3)  to  (x,  y,  z)  it  is  seen  from  eqs.  (2.9) 
—  (2.11)  that  the  x-motion  is  identical  for  all  Ht’ s  and  so  is  now  separable  from  that  of 
y  and  z,  an  anticipated  result  since  the  coordinate  x  was  defined  so  as  to  be 
perpendicular  to  the  plane  containing  all  the  geometrical  changes.  The  calculation  of 
the  dynamics  of  the  electron  transfer  starting  from  electronic  configuration  1,  i.e., 
from  D*BA,  now  involves  only  the  two  coordinates  y  and  z. 

While  eqs.  (2.10)  and  (2.11)  follow  from  (2.5)  and  the  coordinate  transformation 
(2.8),  the  equilibrium  values  in  the  (x,  y,  z)  system  for  the  points  ax ,  a2,  a3 ,  evident 
from  eqs.  (2.10M2.11),  are  also  easily  inferred  from  Fig.  I.7 

III.  The  Rate  Constant 

For  the  loss  of  the  electron  from  D*  we  use  a  nonadiabatic  description  (weak 
interaction  of  D*  with  BA).  A  Golden  Rule  approximation  can  then  be  used  for  a  first- 
order  rate  constant  kr  for  the  disappearance  of  electronic  configuration  1.  For 
reactant  from  a  particular  vibrational  state  of  D*  this  kr  is  given  by3,9 

Ar  =  2n \($N  (y,zt$s  (y,  z))2 -EN  )  ,  (3.1) 

1  23  23  1 

23 

using  the  separability  of  the  r  from  the  (y,  z)  motion.  In  eq.  (3.1),  Hn  is  the 
dimensionless  Hl2  in  (2.4),  denotes  the  energy,  divided  by  hw,  for  a  system  whose 
electronic  configuration  is  1  and  whose  vibrational  wavefunction  is  tyN^y,  2)-  The 
vibrational  wavefunction  for  the  motion  on  a  potential  energy  surface  based  on  the 
pair  of  electronic  configurations  (2,3)  is  denoted  by  $^v23(y,  z),  and  the  corresponding 
vibrational  energy,  divided  by  ft  to,  is  En23  • 

The  stationary-state  Schrodinger  equation  for  ,  obtained  from  (2.9)  -(2.11). 
is 
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(3.2) 


a2  a2 

(-*—  -  +  *(y  +  a1A1/A2)2+  *z2  ](J)^  (y,*)  =  £N  <J>V  (y,z)  , 

dy  dz  111 

Equation  (3.2)  is  separable,  and  we  have 

=  ^  ,  (33) 

where  n  and  m  are  the  respective  quantum  numbers  for  the  y-  and  z-vibrational 
motions  in  electronic  configuration  1,  thus  denoting  the  pair  (n,  m).  4>a(y)  and 
<$>m(z)  are  the  corresponding  wavefunctions. 

When  the  eigenfunctions  $iy23(y, z)  in  eq.  (3-D  approximated  by  being  taken 
to  be  the  solutions  of  the  diabatic  equation, 

/  a2  a2  \ 

(-*71  “*75  +ff22^  +  H22(z)j<l>Vr,^2)  =  EV  ^'2)  ’  0.4) 

dy  dz  mJ  Zo  Zo 


<t>/V23(y.  z)  becomes 

■  V,  =  E .■  +  E„-  ■ 

Z3  Zo 

(3.5) 

where 

£ 

(-*—  +i(y-A1A/a1)2I<|>)|.,2)(y)  -  Hy^  (y)  =  E n,  *  J2' (y) 

(3  6) 

and 

a2 

-i-+H„(z)  <t>m.U)  =  En.*Jz)  , 

(3.7) 

where  Hn(z)  is  defined  in  (2.11).  The  first  half  of  eq.  (3.6)  defines  Hy. 

To  obtain  a  thermally  averaged  rate  constant,  eq.  (3.1)  is  multiplied  by  the 
Boltzmann  factor  exp( -p'Ewjl/Qxy  and  summed  over  Nv  Here,  (!’  denotes  h(x>/kBT,  k3 
being  Boltzmann’s  constant  and  T  the  temperature;  and  Qyt  is  the  yz-vibrational 
partition  function  when  the  electronic  configuration  is  that  of  1.  Since  the  y  and  z 
vibrations  for  the  initial  electronic  configuration  1  are  treated  as  harmonic,  Qsz  is 
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given  by  a  standard  expression.  When  the  y-  motion  is  that  of  a  harmonic  oscillator 
in  the  diabatic  electronic  states  1  and  2,  and  when  the  harmonic  H22(z)  potential  is 
used  instead  of  an  adiabatic  potential  E _{z)  discussed  later,  eq.  (3.1)  yields  the 
standard  result  for  the  thermally-averaged  rate  constant  kr(T)10 

k(T)  =  < kr )  =  2na)|/f12l2(exp(-\12coth(0V2)  0>/2)  )yA12cosech  (072)1  ,  (3-8) 

where  Ip{u)  is  a  modified  Bessel  function  of  order  p  and  argument  u,  and  where 

p=-A£12,  \l2  =  i(Oj  +  a2)  .  (3  9) 

\12  is  the  usual  "reorganization  term  in  electron  transfers,  divided  by  ha,  namely  in 
the  present  case  the  value  of  H22(z)—Hn(z)  —  &El2,  at  the  equilibrium  z-position  for 
the  Hn(z)  potential.11  All  the  quantities  in  (3.1)  and  (3.9),  apart  from  co,  are 
dimensionless. 

In  arriving  at  eqs.  (3.8M3.9)  from  (3.1),  using  a  set  of  z-mode  acceptor  states 
(Hamiltonian  H22(z)),  a  coupling  of  this  mode  to  the  many  other  coordinates  in  the 
system  was  assumed.  Each  z-mode  acceptor  state  was  thereby  broadened.  The  sum 
over  N22  in  (3.1)  can  then  be  replaced  by  an  integral  containing  the  density  (Aco) ~ 1  of 
z-mode  states.  The  condition  Enx  -£V23  =  0  in  the  delta  function  would  be  satisfied 
by  having  a  maximal  number  of  vibrational  quanta  pha  going  into  the  z-mode,  with 
the  residual  energy  imbalance  being  supplied  from  a  coupling  to  a  states.  Thereby, 
the  p  in  (3.9)  is  the  integer  nearest  —  Al?12  (with  -  A£12  expressed  in  units  of  ha).  In 
practice,  eqs.  (3.8)-(3.9)  will  be  treated  as  being  interpolative,  i.  e.,  with  p  not  being 
restricted  to  be  an  integer,  but  rather  equal  to  -  A El2.  The  agreement  of  the  various 
results  in  Part  II,  e.  g.,  a  comparison  of  (3.8)  at  0°  K  with  a  numerical  solution  for 
provides  some  support  for  this  approach. 


IV.  Population  of  B~ 


When  B-  is  formed  as  an  intermediate  in  an  incoherent  (hopping)  mechanism 
its  time-dependent  population  can  be  calculated  by  standard  chemical  kinetics  in 
terms  of  rate  constants  kx  and  k2  of  the  two  successive  steps.  However,  when 
electronic  configurations  2  and  3  are  strongly  coupled  electronically,  the  wavepa^ket 
representing  the  system  can,  on  reaching  the  B-  region  of  the  (2,  3)  pair  of  surfaces, 
move  rapidly  downhill.  It  is  the  dynamics  on  this  (2, 3)  pair  of  surfaces  that  is  treated 
next,  the  packet  being  steadily  resupplied  by  the  nonadiabatic  transition  from  the 
decaying  electronic  configuration  1. 

At  zero  time  the  system  is  in  some  zeroth-order  vibrational  eigenstate  <t>nm(y, 
=  <J>m(z)  of  electronic  configuration  1.  Using  the  type  of  perturbation  theory 

associated  with  exponential  decay  phenomena,8  the  subsequent  time-evolution  of 
$nm(y,  z,  t )  in  electronic  configuration  1  is  given  by3 

-*  t/2-iE  t 

r  nm  4>  (y,z)  .  (4.1) 

nm  firn 

where  the  rate  constant  kr  is  given  by  eq.  (3.1)  and  Enm  is  the  vibrational  energy  for 
the  (n,  m)  vibrational  state  in  electronic  configuration  1. 

The  wavefunction  from  electronic  state  1  is  continuously  deposited  on  the  (2,3) 
pair  of  surfaces  as  a  result  of  the  nonadiabatic  transition.  In  the  perturbation  theory8 
used  to  obtain  eqs.  (3.1)  and  (4.1),  the  deposited  wavefunction  is  |qr2  )H  124>nm(y.  z.  t). 
The  latter  then  moves  in  time  on  the  coupled  (2,  3)  pair  of  surfaces,  evolving  under 
the  influence  of  the  electronic-nuclear  Hamiltonian  H  given  below.  The  electronic- 
nuclear  wavefunction  qr(y,  z,  t)  is  now  given  by  the  convolution  expression,12 

f*  -iH  U-9) 

lv(y,z,  0)  =  -*  j  dQe  n  \y2)  Hl2$nm(y,z,Q)  .  (4  2) 

where  <t>nm(y,  z,  0)  is  given  by  eq.  (4.1).  It  is  assumed  for  simplicity  in  eqs.  (4.2)  and 
(3.1)  and  that  because  of  the  large  separation  distance  the  electronic  matrix  element 
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Hl3  directly  coupling  state  qr,  to  qr3  can,  in  the  first  approximation,  be  neglected.  In 
terms  of  a  bra-ket  notation  for  the  diabatic  electronic  wavefunctions  qr’s,  the  H  in 
eq.  (4.2)  can  be  written  as 


H 

y* 


+  H 

i 


(4.3) 


where,  without  loss  of  generality,  the  definition  of  Hz  includes  the  H23  term: 


(4.4) 


and,  using  the  fact  from  eq.  (2.11)  that  H22(y)  =  if33(y), 

»,  =  i  i  • 


(4.5) 


3 

Here,  /  denotes  the  identity  operator  2  I'PiX  'Pj  in  the  (xp2,  qr3)  subspace,  and 

i  =  2 

Hy'  is  defined  in  the  first  half  of  (3.6):  Hy  is  seen  from  eq.  (4.5)  to  induce  no  transition  s 
in  the  (qr2,  qr3)  subspace,  as  expected. 

Since  Hy  and  Hz  commute,  the  time-evolution  operator  exp  (  —  iHyzx )  in  eq.  (4.2) 
can  be  written  as  exp  ( —  iHyx)  exp  ( —  iHzx).  Further,  using  (4.5)  we  have 

exp  ( — i^f^x)  |tp2  ><J>fiCy)  =  |\p2>exp(-iH^)<|>ri(y)  ^ ^ 


With  a  change  of  variable  f  —  0-»0,  eq.  (4.2)  now  becomes 


-(*  f2+iE  )t  (*  /2+t E  )0  -i H  0  - iH  0 

l'p(y.*,<)>=  r  nm  dQe  r  nm  (e  y  4>a(y)He  *  >  (4 .7) 


which  is  the  desired  final  expression. 

The  "population”  B_(0  of  the  intermediate  electronic  configuration  2  is 
obtained  by  projecting  |q/(y,  z,  *)>  onto  the  diabatic  electronic  state  |ip2  )  and 
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integrating  over  the  nuclear  coordinates: 


B'tt) 


■ 

dy  dz\{y2 1  z.rt)  |2 


(4.8) 


Strictly  speaking,  B  ~  does  not  actually  exist  in  this  mechanism,  and  it  would  be  more 
rigorous  to  treat  an  observable,  such  as  the  absorption  spectrum  in  the  B  spectral 
region,  rather  than  defining  aB'(t)  via  eq.  (4.8).  Use  of  eq.  (4.8)  for  B-  presumes  that 
the  quantity  so  defined  adequately  desciibes  this  depletion  in  the  spectral  region  of  B. 

The  exp(  —  iH  '9)$„(y)  term  in  eq.  (4.7)  can  be  written  as 

e"  y  ♦„<*>=  Z  <$„.<2>  &)!♦„<?»«  “*  4>n.(y)  .  (4.9) 


where  4>„.(2)(y)  is  the  n’  th  eigenfunction  of  the  Hamiltonian  H '  defined  in  eq.  (3.6), 
En,  is  the  corresponding  eigenvalue,  and  <$,/2)(y)|<l>„(y))  is  a  Franck-Condon  factor. 
Each  of  the  quantities  in  eq.  (4.9)  is  well-known13  for  the  present  case  where  Hy  is  a 
harmonic  oscillator  Hamiltonian,  eq.  (3.6).  The  z-dependent  term  in  (4.7),  exp  ( - 
i/fz0)|qr2>  4>^(2),  can  be  calculated  using  a  Fast  Fourier  Transform  (FFT)  method.3  In 
the  procedure  some  mechanism  is  prescribed  for  the  dissipation  of  the  wavepacket 
after  reaching  the  H33(z)  surface,  for  example  by  adding  an  imaginary  term  —  ir  to 
H33(z)  in  eq.  (4.4)  or  by  modifying  H33(z)  and  then  introducing  an  absorbing 
boundary.4 

In  ref.  4,  for  comparison  with  some  of  the  numerical  results  for  eqs.  (4.7M4.8),  we 
have  also  made  some  calculations  with  H22iy)  +  H22  (z)  in  eq.  (3.4)  replaced  by  an 
adiabatic  surface  E_(y,  z).  From  a  comparison  of  the  two  sets  of  calculations 
information  is  obtained  on  the  reflection,  if  any,  of  the  wavepacket  during  its  passage 
through  the  H22  and  H33  intersection  region.  E _{y,  z)  is  the  lower  of  the  two  adiabatic 
potentials  constructed  from  H22  and  H33 ,  and  is  given  by 
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E_(y,z)  =  HHg+^-[(H2-^2+4^ 


(4.10) 


where  HJ*  is  Hu(y)  +Hu(z),  defined  in  (2.10M2.11). 

Inasmuch  as  H22y*—H3J*  depends  only  onz(cfeq.  (2.11))  the  adiabatic  potential 
EJy,  z)  given  by  (4.10)  is  seen  to  be  the  sum  of  two  terms, 


£_(y,z)  =  *(y-AlA2(a1)2  +  EJz) 


(4.11) 


where 


EJz)  =  *  {[H^z)  +  H 33(2)] -  H^Jz)) 2  +  4H^}* }  (4.12) 


V.  Energy  Partitioning  Among  the  y-  and  z-  modes 

When  the  wavepacket  from  state  1,  given  in  (4.2),  is  deposited  on  the  (2,  3) 
surface,  the  population  B-(0  estimated  in  (4.8)  depends  on  the  partitioning  of  the 
excess  energy  among  the  y  and  z  modes.  To  obtain  some  insight  into  numerical 
results4  for  this  partitioning,  and  later  into  an  effective  rate  constant  kv  it  is  useful  to 
introduce  here  a  quite  rough  but  analytic  estimate: 

The  intersection  of  the  Hn(y,  z)  and  H22(y,  z)  potential  energy  surfaces  has  some 
lowest  point  (y12*,  zl2%).  Using  the  arguments  given  in  Appendix  A,  the  excess  energy 
Ey  in  the  y-mode,  measured  relative  to  the  energy  at  the  lowest  point  y2  of  the  curve 
tf22(y),is 


Ey  =  n  +  t+  (A £12/X12)(X2/X23)- A£12-n*  +  i  (5.1) 

where  this  n'  so  defined  is  predicted  to  be  a  dominant  term  in  eq.  (4.9)  or  an  average, 

E  n'l<<b  ,<2>l<t>n)l2,  of  the  terms  there.  Subtracting  this  E  from  the  approximate  total 

n'  y 

energy  after  the  ipj-npj  transition,  m+n  + 1  —  Ai?12,  the  average  z-mode  energy  Ez  for 
motion  on  the  (2, 3)  pair  of  surfaces  is  estimated  to  be 


m  +  i-WVA12)(X,2%)  ■  <5  2> 

In  the  three-state  dynamics  the  system  is  initially  deposited  from  state  1  onto  the 
H22(y,  z)  surface,  since  Hl3  —  0.  However,  because  of  the  strong  coupling  of  the  (2,  3) 
pair  it  may  quickly  settle  down  to  a  motion  largely  on  the  E_(y,  z)  surface. 
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In  that  case  there  is  a  larger  energy  available  at  (y2°,  z2°),  largerer  than  (5.1)- 
(5.2)  and  distributed  between  the  y-  and  z-  motions.  The  additional  amount  of  energy 
shared  between  those  coordinates  is  estimated  in  appendix  A  to  be  H22(z2°)-E  _(z2°), 
and  thereby  that  Ey  +  Ez  is  given  by 

£  +  £  =  m  +  n  +  1  -  &El2  +  ti(\n  +  iEnf  +  (5.3) 

instead  of  m  +  n  + 1  —  AE12. 

To  interpret  some  of  the  numerical  results  for  B~(f)  and  for  given  in  Part  II, 
and  provide  an  approximate  indication  of  what  might  be  expected  for  other  values  of 
the  various  parameters,  we  obtain  below  a  rather  rough  classical  estimate  of  the  time 
spent  by  the  system  in  the  "B~  region”  and  for  an  effective  rate  constant  k2 .  For  this 
purpose  we  approximate  the  length  of  this  region  by  the  z-interval  between  z12*  and 
the  intersection  z23*  of  the  ff22(z)  and  ff33(z)  surfaces.  The  point  z23*  is  given  by 
equations  analogous  to  eqs.  (A3)  and  (A4): 

4>=2;+"(2;-2S>  •  (54> 


where  M  satisfies 


-(2M+l)A2J  =  i£23  (5  5) 

Thereby,  the  value  of  z2, 3*  -  z12* ,  obtained  as  the  difference,  z23*  —  z2°  —  (z12*  -  z2°),  is 
found,  using  eqs.  (A3)  and  (A4)  for  z12*,  to  be 

4  -  4  -  +  V  +  -  »nXV‘u»  '« V*  (S  6> 

When  the  wavepacket  undergoes  negligible  reflection  in  the  region  near  z23\  a 
condition  we  denote  by  setting  the  2-»3  transition  probability,  w23,  to  be 
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approximately  unity,  the  time  t2  spent  by  the  system  in  the  B-  region  is  then, 
classically, 

W."?  (when  ^23  ~1)  (57) 

*12  * 

where  vz  is  the  local  2-component  of  velocity.  In  the  case  of  the  H22(z)  curve  vz  is 
[2 {Ez—H22(z)+H22(z°2)}]*,  and  in  the  case  of  the  E  _(z)  curve  vz  would  be  somewhat 
larger.  Eq.  (5.7)  can  be  written  in  terms  of  some  averaged  velocity  0z ,  t2  =  (z*23- 
z*12)/iV 

If  t2  is  regarded  as  the  reciprocal  of  an  effective  rate  constant  k2,  then  eq.  (1.1)  still 
applies  approximately  for  the  maximum  population  of  B~,  but  now  with  the  k2  given 
by 


k2~uz' l*a  -  •*«!  (when  —  1)  (5-8) 

A  few  examples  of  B-^  calculated  from  eqs.  (1.1)  and  (5.8)  are  given  in  Part  II,4  and 
indicate  that  eq.  (5.8)  is  correct  to  roughly  a  factor  of  about  two  when  u/23  is 
approximately  unity.  Further,  a  comparison  is  also  made  there  with  the  full  B~U)  vs. 
t  curve  using  this  effective  k2  .  This  k2  becomes  large  when  the  z-interval  in  the 
denominator  of  eq.  (5.8)  becomes  small  (but  not  too  small,  if  the  expression  is  to  have 
some  validity),  and  if  at  the  same  time  u/23~l. 

In  several  respects  a  recent  article16  by  Lin  on  a  donor-bridge-acceptor  electron 
transfer  is  complementary  to  the  present  one:  the  case  when  the  B~  orbital  is  close  to 
D*  is  included  there  and  the  overall  D*-»A  transfer  is  assumed  to  be  coherent.  There 
are  several  differences:  a  steady-state  (Golden  Rule)  perturbative  type  treatment 
was  used  in  ref.  16  for  the  entire  process,  D*-*A,  so  that  the  time-evolution  for  the 
formation  of  A-  was  described  by  a  single  exponential.  The  electronic  matrix 
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elements  are  presumably  assumed  to  be  small  in  ref.  16,  so  that  the  perturbative 
treatment  can  be  used  throughout.  In  the  present  treatment  a  highly  nonexponential 
formation  of  A~  may  occur  (e.g.,  Fig.  3  in  Part  II),  and  a  large  ff23  matrix  element  is 
permitted.  One  focus  in  the  present  article,  unlike  that  in  ref.  16,  is  on  the  maximum 
"B~”  population.  Another  treatment  akin  to  that  in  ref.  16  is  that  of  Kharkats  et 
a/.,17  who  used  for  simplicity  a  one-coordinate  treatment.  Other  studies  include  one 
by  Friesner  and  Wertheimer,18  who  assumed  a  strong  vibronic  (nonperturbauve) 
coupling  in  a  photosynthetic  system:  They  employed  vibrational  frequencies  which 
nearly  matched  the  assumed  electronic  energy  gaps  and  so  caused,  thereby,  a  strong 
electronic  interstate  mixing. 

Applications  of  the  present  formalism  are  made  in  Part  II. 
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Appendix  A.  Energy  Partitioning  Among  the  x-  and  y-  Modes 

The  average  distribution  of  the  excess  energy  in  each  of  the  y  and  z  modes  after  a  1 
-»  (2,3)  transition  is  estimated  in  this  Appendix.  When  a  matrix  element  (<t>N{  (y, 
z)\<t>N23(y,  z))  in  eq.  (3.1)  is  treated  semiclassically,14  it  has  a  stationary  phase  point 
which  lies  on  the  intersection  surface  of  Hn(y,  z)  and  a  second  surface,  H22(y,  z)  or  E 
(y,  z),  depending  on  the  approximate  $jv23(y,  z)  used  for  the  stationary  phase 
condition:  The  stationary  phase  point  occurs  where  there  is  no  change  in  the  y-  and  z- 
components  of  momentum  and  so  occurs  on  the  above  intersection  surface  (Franck- 
Condon  principle).  In  the  classical  limit,  the  dominant  semiclassical  matrix  elements 
provide  information  on  an  ensemble  of  classical  trajectories  from  state  1  to  state  (2,3). 
This  ensemble,  in  their  crossing  from  the  Hn(y,  z)  surface  to  the  second  surface,  is 
centered  around  the  lowest  point  ( y12* ,  zl2)  on  the  intersection  surface. 

If  the  y-mode  and  2-mode  energies  before  and  just  after  the  transition  are  En  and 
Em,  respectively, .in  the  notation  of  eq.  (3.3),  the  y-mode  kinetic  energy  at  (yl2*,  z12*)  is 
En  -  i  (y12*  -y°),2  where  (y^,  Zj°)  specifies  the  minimum  of  the  ifn(y,  z)  surface.  The 
potential  energy  for  the  y-mode  on  the  H22  (y,z)  surface,  relative  to  what  it  would  be  at 
the  point  y2°  is  i(y12*  -y20),2  where  (y2°,  z2°)  specifies  the  minimum  of  H22  (y,  z).  Thus, 
Ey  (  the  y-mode  energy  in  excess  of  the  potential  energy  if22(y2°)  after  the  1  -*  (2,3) 
transition,  is 

£,  =  »  +  i-*(y{2-yp2  +  i(y\2-yp2  •  (A1) 

We  define  the  corresponding  z-mode  energy  Et  as  its  value  relative  to  that  at  z2°. 
Thereby, 

Ez  =  m  +  i-i(z\2-e^2+  *  (z\2  -  zf  (ffa  surface)  (A2) 
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The  dominant  n'  term  in  eq.  (4.9)  is  the  one  whose  energy,  n'  +  in  units  of  Aco,  most 
closely  approximates  the  energy  Ey  given  by  (Al). 

Using  concepts  employed  for  intersecting  parabolic  surfaces15  the  point  (y12*,  z12*) 
lying  on  the  intersection  of  Hniy,  z)  and  H22{y,  z )  surfaces  and  having  the  least 
potential  energy  is  given  by 

y\2=y°i  +M(y°i  •  4=z?+M(z?-^)  .  (A3) 

where  Af  satisfies  the  equation15 

-(2Af+  1)X12  =  A£12  .  (A4) 

Equations  (A1)-(A5)  then  yield  the  energy  partitioning  expressions,  eqs.  (5.1M5.2). 

If  the  system  after  the  l-»3  transition  is  better  described  near  z2°  as  residing  on 
the  E_{y,.z)  surface  rather  than  on  the  H22(y^z)  one,  then  there  is  available  for 
distribution  among  the  y-  and  z-  motions  some  extra  energy  (extra  when  H22>E _), 
namely,  H22{z  2°)  —  E_  ( z  20),.  Using  eq.  (4.12)  and  (2.11)  this  additional  energy  term 
is  found  to  be  given  by  eq.  (5.3). 
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Appendix  B.  Landau-Zener  Transmission  Probability  Near  z23*. 

In  a  Landau-Zener  treatment  the  probability  w23  of  remaining  on  the  potential 
energy  curve  E_(z )  after  a  single  passage  through  z23t  is  given  by14 

-2nH2/u\tU\ 

Til  -  1  _  CX  Z  _  v 


when  z23t  is  real.  Here,  vt  is  the  velocity  at  the  crossing-point  and  As  is  the  difference 
of  slopes  of  the  intersecting  diabatic  potential  energy  curves  there.  If  we  set  vz 2  =  N 
in  dimensionless  units,  (which  defines  N),  and  if  we  introduce  the  value1  for  |As|, 
namely  (z3°  —  z2°)  in  dimensionless  units,  one  finds  that 


wn=  1 


-2n"*2/SX'2 


(B2) 


For  example,  if  the  z-energy  in  the  vicinity  of  z23*  is  approximately  the  zero-point 
energy,  then  1.  Employing  eq.  (B2)  when  H23  becomes  very  large  may,  however, 
be  questionable. 
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Legend  to  Figure 

Fig.  1  Original  (q{  q2  q3)  and  rotated  ( x  y  z)  coordinate  axes.  The  x-axis  is 
perpendicular  to  the  plane  containing  at  a2a3  and  intersects  with  it  at  Q.  The 
z-coordinate  axis  is  parallel  to  the  line  a2  a3  and  serves  as  the  reaction 
coordinate  after  the  initial  loss  of  the  electron  from  D*.  The  y'  axis  is  the  result 
for  the  new  g2-axis,  obtained  after  the  rotation  of  an  angle  $  about  the  q^-axis, 
and  the  y-axis  results  from  Oy'  after  the  rotation  about  0 z  through  an  angle  0. 
P  is  the  foot  of  the  perpendicular  from  0  to  the  line  a2  a3. 


